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^Netherlands Consortium for Systems Biology, Netherlands Institute for Systems Biology, Amsterdam, 
The Netherlands 
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A key step in blood vessel development (angiogenesis) is lumen formation: the 
hollowing of vessels for blood perfusion. Two alternative lumen formation 
mechanisms are suggested to function in different types of blood vessels. The 
vacuolation mechanism is suggested for lumen formation in small vessels by 
coalescence of intracellular vacuoles, a view that was extended to extracellular 
lumen formation by exocytosis of vacuoles. The cell- cell repulsion mechanism is 
suggested to initiate extracellular lumen formation in large vessels by active 
repulsion of adjacent cells, and active cell shape changes extend the lumen. 
We used an agent-based computer model, based on the cellular Potts model, 
to compare and study both mechanisms separately and combined. An exten- 
sive sensitivity analysis shows that each of the mechanisms on its own can 
produce lumens in a narrow region of parameter space. However, combining 
both mechanisms makes lumen formation much more robust to the values of 
the parameters, suggesting that the mechanisms may work synergistically 
and operate in parallel, rather than in different vessel types. 



1. Introduction 

Blood vessels are essential for efficient transport of oxygen and nutrients through- 
out the body. During embryogenesis, a system of blood vessels is formed. 
Remodelling of the vasculature continues throughout our entire life. Sprouts of 
endothelial cells grow out from the existing vasculature to reach oxygen-deprived 
regions, including wounds or tumours, in a process called angiogenesis [1]. Sub- 
sequently, the new sprout hollows to allow blood to perfuse [2,3]. This hollowing 
is called lumen formation, and can occur in the absence of blood pressure, such 
as in vitro [3-6] and also in vivo in intersegmental vessels (ISVs) of zebrafish [6]. 
Which mechanisms are responsible for lumen formation is debated [7-12]. 

Lumen formation is extensively studied in epithelial tissues, which has 
resulted in a range of potential mechanisms [13]. Three of these can form 
lumens in cords of cells: cavitation, cell hollowing and cord hollowing. Cavita- 
tion is unlikely to play a role in lumen formation of endothelial tubes, where 
apoptosis is rarely observed [3]. The two remaining mechanisms assume 
intracellular lumen formation within cells in unicellular tubes (cell hollowing) 
versus extracellular lumen formation between cells in multicellular tubes 
(cord hollowing). The debate whether lumens form intracellularly [14] or extra- 
cellularly [15] in blood vessels originates from the nineteenth century [16]. This 
led to two opposing views on the molecular mechanisms of lumen formation in 
endothelium: the vacuolation mechanism [3,4] and the cell -cell repulsion 
mechanism [5,17]. The vacuolation mechanism (figure \a,h) originally suggested 
that an intracellular lumen is formed by coalescence of vacuoles [3,18], and 
this view was extended with extracellular lumen formation by exocytosis of 
vacuoles [6,19]. The cell-cell repulsion mechanism (figure Ic) assumes that extra- 
cellular lumens initiate by active repulsion of adjacent cells and are expanded 
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Figure 1. Experimental observations of lumen formation, (a) Time-lapse images of endothelial cells (ECs) in a three-dimensional collagen gel, showing formation of 
vacuoles into larger intracellular compartments. (Adapted from Kamei et al. [3].) (b) Two-photon time-lapse imaging of an ISV in which membranes are labelled 
with EGFP-cdc42wt and intravascularly injected red quantum dots are serially transferred by vacuole fusion. (Adapted from Kamei et al. [3].) (c) Lumens form during 
in vitro three-dimensional angiogenic sprouting assays by cell -cell repulsion facilitated by negatively charged CD34-sialomucins (control panel). Cleavage 
(neuraminidase panel) or neutralization (protamine sulfate panel) hereof reduces lumen formation. (Adapted from Strilic et al. [5].) 



by active cell shape changes [5,17]. As vacuoles are often 
observed in ISVs, but not in aortae, the vacuolation mechan- 
ism is suggested to function in small vessels and cell -cell 
repulsion is suggested to function in large vessels [6,9-12]. 

It is difficult to distinguish the two proposed lumen for- 
mation mechanisms experimentally, because they use similar 
proteins and pathways — everything is intertwined. Therefore, 
we used a computational model, which allows us to isolate, 
modify and study single mechanisms and selected components 
[20-22], and compared the efficiency of both lumen formation 
mechanisms separately and combined. An extensive parameter 
sensitivity analysis of the model suggests that lumen formation 
is most robust to inhibitions of either mechanism when the two 
mechanisms are combined, indicating that they function syner- 
gistically in lumen formation. Thus, our model challenges the 
view that each of the mechanisms operates on its own in differ- 
ent types of blood vessels [6,9-12], and supports the idea that 
the different mechanisms operate in parallel [7]. 



2. Results 

We developed a two-dimensional, multi-scale, agent-based 
computational model in which lumen formation emerges 
from predefined behaviour of components at the cellular and 
subcellular scale. During angiogenesis, lumens form shortly 
after new sprouts have originated [2,3]. Thus we can assume 
that sprouting and lumen formation are separated in time, 
such that we can model lumen formation in a preformed 



sprout. The sprout is represented as a branched cord of 12 
cells within an extracellular matrix (ECM; figure 2a, b). The 
bifurcating geometry contains narrow, single-cell-wide regions 
at the tips of the branch and two-cell-wide regions at the trunk 
of the branch. The model currently neglects degradation or 
secretion of ECM by the cells that we considered in previous 
work [23]. ECM fluid is present at each end of the branch to 
allow stretching rather than widening (as seen in experiments; 
figure la,c) of the vessel during lumen formation. 

The model represents a two-dimensional cross section of the 
vessel. In a cross section, a cell with an intracellular lumen will 
be split up, whereas it would be doughnut-shaped in three 
dimensions. Despite this obvious disadvantage of using a 
two-dimensional simulation, it still gives us useful insights 
into the local mechanisms driving lumen formation. Further- 
more, the computational efficiency of the two-dimensional 
model allows us to evaluate the behaviour of the model for a 
large number of parameter settings. 

The model consists of four modules (figure 3) that we 
combine differently to represent each of the mechanisms. Cell 
motility is regulated by the cellular Potts model (CPM) [1A,15\, 
which considers the shape of cells and their adhesive properties. 
Agent-based extensions of the CPM represent subcellular struc- 
tures, including membrane proteins, vesicles and vacuoles. Cell 
surface 'polarization results in a basolateral membrane that lines 
the vessel and connects the ceUs, and an apical membrane 
where the lumen will form. During vacuolation, pinocytotic ves- 
icles are formed within the cells, coalesce into vacuoles and are 
secreted at the apical membrane. Cell- cell repulsion occurs by 




Figure 2. Model initialization, (a) The configuration of the model at initialization, (b) Final outcome of a simulation without cell surface polarization, (c) Final 
outcome of a simulation with cell surface polarization. 
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Figure 3. Modules of the computational model, {a) Cell motility by a copy of the membrane outwards (1) and inwards (2). (b) Cell surface polarization, resulting in 
cytoplasm (blue) surrounded by a basolateral membrane (grey) and an apical membrane (orange), (c) Vesicles (red) can be created by pinocytosis (3), move in the 
cell by a biased random walk (4), fuse into vacuoles (yellow) (5) and be secreted (6) into luminal fluid (light blue), (d) The apical membranes of different cells 
repulse each other (7) owing to the negative charges on the extracellular domains of the CD34-sialomucins herein. 




active repulsion of apical membranes from opposing cells. We 
briefly describe each module here; for details and reference 
settings (tables 1-3) see Material and methods. 

2.1. Cell motility module 

The CPM [24,25] is a convenient model to deal w^ith structural 
and spatial aspects in lumen formation, because it considers 
the shape of cells and their physical interactions w^ith their 
surroundings. Cells are represented as patches of connected 
sites of a lattice. The diameter of a cell is approximately 
5-10 |jim (one pixel = 250 nm), and a cell is initialized w^ith 
625 pixels ( + 25 x 25). We represent the subcellular scale 
using a compartmentalized CPM [26-28]. Each lattice site x 
in a compartment is associated w^ith the same unique compart- 
ment identifier {cr{x)). Depending on the mechanism that is 
tested, compartments w^ith the type {T{a{x))) cytoplasm, apical 
membrane, basolateral membrane, vesicle and vacuole can exist 
w^ithin a cell. ECM, ECM fluid and ceUs have an additional 
unique cell identifier ^{a{x)), w^ith the set of compartments a 
that belong to the same cell having the same cell identifier. 
The ECM is immobile, but cells and compartments move 
by copying pixels at the membrane inw^ards or outw^ards 
(figure 3a). The source pixel (x) together w^ith a random, adja- 
cent target pixel {x') is randomly selected for each copy attempt, 
resulting in a stochastic simulation. Acceptance of these pixel 
copies depends on properties of the cells and compartments, 
such as adhesion/ repulsion strength (contact) and size (area). 
Adhesion or repulsion strength can result from protein inter- 
actions and from surface tensions on membranes and is 
modelled by contact energy. A distinction is made betw^een 



internal contact energy betw^een compartments of the same 
cell (/i(t,t)) and external contact energy between compartments 
of different cells (/e(t,t)). As a reference contact energy, we use 
a value of / = 10: a low^er contact energy leads to adhesion, 
whereas a higher contact energy leads to repulsion. Lattice 
sites that are not occupied by ECM or cells are of type ECM 
fluid or luminal fluid. ECM fluid is already present at initializa- 
tion next to the branch and is in contact with the ECM. Luminal 
fluid is secreted by vesicles and vacuoles in between cells. 
Luminal fluid will fuse to the ECM fluid when surrounded 
by ECM fluid or when in contact with the ECM. Cells, vesicles, 
vacuoles and subvolumes of fluid have a preferred size (A), 
which is conserved on average, with the elasticity parameter 
>^type (see Material and methods for details) regulating the 
allowed deviation from the preferred size [24]. For cells, ves- 
icles and vacuoles, this reflects the semipermeability of 
membranes for water. For subvolumes of fluids, this results 
in near incompressibility, which could resemble a hydrostatic 
pressure in the sublumens. 

2.2. Cell surface polarization module 

Cell surface polarization into a basolateral and an apical 
domain is the first step in lumen formation and is regulated 
by integrin signalling from the ECM [9,17,29-33]. In the com- 
putational model, the membrane (re)polarizes every other 
time step into an apical and a basolateral compartment 
based on the relative position to the ECM (figures 2c 
and 3b). During repolarization, each 'mispolarized' pixel 
in the membrane is assigned to the correct compartment. 
A membrane pixel becomes part of the basolateral 
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Table 1. Reference values of external contact energy. The external contact energy (J^) is listed for each type combination. 
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Table 2. Reference values of internal contact energy. The internal contact energy (7|) is listed for each type combination. 
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Table 3. Reference parameter values. The reference value is listed for parameters that relate to probabilities, the motility /jl and the elasticity A for certain types. 

/^A=1.0 /'fuse =1.0 /'pin = 1.0 /X = 50 Ace|| = 7 Afluids = 6 Avacuole = 50 Avesicle = 1000 



compartment when it is in direct contact with the ECM. To 
form the lateral junctional regions of the cells, membrane 
pixels that have at least two neighbouring pixels that are in 
direct contact with the ECM also become part of the basolat- 
eral compartment. The rest of the membrane pixels become 
part of the apical compartment. Contact energy is used to 
mimic adherens junctions (Ji {apical, basolateral)) and to set 
the surface tension of the apical membrane (/e (fluids, apical)). 



2.3. Vacuolation module 

The vacuolation mechanism is one of the two proposed mech- 
anisms for lumen formation in blood vessels [3,18]. Vacuoles 
are often observed in vitro and in vivo during lumen formation 
as summarized by Davis & Bay less [4]. Vacuoles were 
observed with electron microscopy (figure la) and by 
expression of the green fluorescent protein (GFP) fusion pro- 
teins GFP-Racl, GFP-Cdc42 and Moesinl-EGFP, which 
colocalize with the vacuoles [4,34-36]. Vacuoles show highly 
dynamic behaviour, as they continuously fuse together, 
grow, shrink and disappear [3,6]. Kamei et al. [3] showed that 
a label (carboxyrhodamine) added to the medium is taken up 
into the vesicles by pinocytosis and is transferred to vacuoles 
by fusion of vesicles into vacuoles. Fusion might be facilitated 
by caveolin-1 because it concentrates at vacuole -vacuole con- 
tact areas [4]. Pinocytotic vesicles are probably trafficked 



along microtubules and actin filaments [37]. Targeting to the 
apical membrane might involve Cdc42 and Moesinl [4,6], 
which both co-localize with vesicles and have a high affinity 
for phospholipids specific for the apical membrane. Even- 
tually, vesicles and vacuoles bridge the entire cell [3] and are 
exported at the apical membrane through exocytosis to create 
a fluid-filled space for a lumen [19]. Each of these steps is 
explicitly modelled in this module. 

To mimic pinocytosis, we assumed that if a membrane 
pixel becomes internalized in the cell owing to cell movements, 
it has a probability (Ppin) to become a new compartment of 
iype vesicle (figure 3c; step 3) or to become part of the cytosol. 
The vesicle diameter is approximated at 250 nm [38], which is 
equal to one pixel length in the model. This size is achieved by 
means of a target area of 1 and a high inelasticity (Aveside)- Ves- 
icle transport is precisely regulated with vesicle-associated 
motor proteins that walk along microtubules and actin fila- 
ments [37], which gives rise to stochastic motion of vesicles 
if the cytoskeleton is randomly oriented [37-39]. Random 
binding and unbinding of vesicles and associated motor pro- 
teins with a randomly oriented cytoskeleton leads to 
diffusive transport behaviour [39]. In our model, we therefore 
model vesicle transport by a random walk with stepping prob- 
ability Pa (see Material and methods), biased by preferential 
adhesion to the apical membrane and fusion into vacuoles 
(figure 3c; step 4). Vesicles preferentially adhere to the apical 



membranes, vesicles and vacuoles by considering their contact 
energy in their stepping. Fusion of neighbouring vesicles, 
vacuoles or a combination of the two happens with probability 
Pfuse (figure 3c, step 5) and generates a single vacuole compart- 
ment for which the target areas are combined. A vacuole moves 
by means of the usual CPM rules (see Material and methods) 
and is only restricted by its contact energy and size constraint. 
When a single pixel of type vacuole is no longer in contact with 
vesicles or vacuoles, it becomes a vesicle. Secretion occurs 
when a vesicle or vacuole pixel is in the apical membrane 
(figure 3c; step 6). This pixel, together with all connected 
pixels of type vesicle or vacuole, is secreted to form a fluid- 
filled luminal space of type luminal fluid. The cell membrane 
is subsequently repolarized, which also leads to an increase 
in the size of the apical membrane. 

2.4. Cell -cell repulsion module 

The alternative cell -cell repulsion mechanism [5,17] assumes 
active, electrostatic repulsion of the apical membranes of adja- 
cent cells, followed by active cell shape changes to extend the 
lumen. Previously, it was suggested that lumens can form by 
relocalization of junctional complexes over the cell membrane, 
leading to local differences in cell adhesion at the future 
luminal and basolateral sides of cells [40-42]. In addition to 
such differential adhesion-driven lumen formation, Strilic 
et al. [5,17] showed in the developing mouse aorta and in a 
three-dimensional angiogenic sprouting assay that apical 
membranes actively repulse each other during lumen for- 
mation by expression of CD34-sialomucin glycoproteins, 
such as PODXL. These transmembrane glycoproteins have 
negatively charged extracellular domains and are transported 
to the apical membrane in preformed vesicles. 

This active repulsion is modelled with high contact energy 
between apical membranes of adjacent cells (/e {apical, apical); 
figure 3d; step 7). Pixels of type cytoplasm can be at the mem- 
brane, because cell surface polarization is only performed 
every other time step in the model to allow some freedom of 
movement for the membranes and for internalization of mem- 
brane pixels for pinocytosis. Therefore, /e {cytoplasm, apical) is 
set to the same value. Je {apical, apical) and Je {cytoplasm, 
apical) together will now be referred to as /rep. Note that this 
module of cell -cell repulsion only includes short-range, 
electrostatic membrane repulsion owing to CD34-sialomucin 
glycoproteins. The additional mechanisms in our model, 
including cell shape changes and invasion of luminal fluid, 
may act to further separate the membranes from each other. 

2.5. Reconstruction of the mechanisms 

Each lumen formation mechanism can be reconstructed 
from the four modules: cell motility, cell surface polarization, 
vacuolation and cell -cell repulsion. The first three modules 
are combined to study the vacuolation mechanism, the first 
two and last for the cell -cell repulsion mechanism and all 
four modules together to study them combined, 'the combined 
mechanism'. A simulation is run for 10 000 time steps to allow 
the creation of continuous lumens. During one time step, the 
Monte Carlo step (MCS), as many copy attempts as there 
are pixels in the grid are attempted. Continuous lumens are 
formed in approximately 24 h [4,5]; a time step in the model 
thus corresponds to a few seconds in real time. 

Because the model has a strong stochastic component, we 
run each simulation 30 times. Efficiency of lumen formation is 



expressed as the fraction of simulations that have formed a con- 
tinuous lumen at the end of a simulation (continuity fraction). 
The lumen is continuous when the fluid in the three tips 
of the branch is part of one single connected component. 
Figure 4 shows a time-lapse image for each mechanism at refer- 
ence settings (corresponding to the electronic supplementary 
material, videos SI -S3). During cell-cell repulsion (e.g. figure 
4a, MCS 10 000), the cells stretch and flatten along the ECM, 
very similar to the experimental pictures in figure Ic. The 
shape and size of the sub lumens that are formed during vacuo- 
lation (e.g. figure ^h, MCS 500) visually resemble the 
experimental observations of Kamei et al. [3] (figure la). 
Although lumens do not form by vacuolation at reference set- 
tings, they do form for higher pinocytosis rates. The model is 
robust to practically all changes in parameters that describe 
the transport and fusion of vesicles and vacuoles (see electronic 
supplementary material, figure SI and text SI). 

Ideally, all parameters in the model would have a quantitat- 
ive value derived from experiments to allow comparison. 
However, most values are not known experimentally. Notably, 
Klann et al. [43] point out that they could not quantify time con- 
stants of reactions and transport processes in vesicle d3mamics, 
because present experimental results focus on the functional 
and qualitative identification of molecular interactions and path- 
ways rather than on the d3mamics of the system. Fortunately, our 
computational model can still be validated in a qualitative 
manner. Most parameters in the lumen formation model corre- 
spond to a protein or molecular process in lumen formation. 
As a result, qualitatively reducing the value of a parameter cor- 
responds to a molecular knockout or inhibition of the molecular 
mechanism it represents and should produce similar effects. 

2.6. Vacuolation requires a high pinocytosis rate 

In experiments, inhibition of pinocytosis prevents lumen 
formation and vacuole formation. Pinocytosis is an integrin- 
dependent process and is inhibited during these experiments 
by blockage of these sites (a2|31 in collagen matrices [18] or 
avp3 and a5|31 in fibrin -fibronectin matrices [44]). 

To mimic this blockage, we reduced the rate of pinocytosis 
(Ppin) in our model. In agreement with the experiments, con- 
tinuous lumens formed by vacuolation for high pinocytosis 
rates, but not for low rates (figure 5a, h). 

2.7. Cell -cell repulsion requires a high cell -cell 
repulsion strength 

Reduction of the cell -cell repulsion strength by cleavage or neu- 
tralization of the negative extracellular charged domains of 
CD34-sialomucins also inhibits lumen formation [5]. Lumen for- 
mation is rescued by subsequent addition of negatively charged 
dextran sulfate that binds to cell surfaces. 

To mimic the cleavage /neutralization in our model, 
we reduced the repulsion strength (/rep). As seen in the exper- 
iment, continuous lumens formed by cell -cell repulsion 
for high repulsion strengths, but not for low values of this 
parameter (figure 5c, e). 

2.8. Vacuolation and cell -cell repulsion combined 
produce lumens most robustly 

Figure 5h,c shows that lumen formation by the combined 
mechanism (green curve) is more robust to changes in the 




pinocytosis rate than the vacuolation mechanism (red curve) 
and that the combined mechanism is more robust to changes 
in the repulsion strength than the cell -cell repulsion mechan- 
ism (blue curve), respectively. The combined mechanism 
already has a continuity fraction of 0.8 in the absence of pino- 
cytosis (figure 5h). The continuity fraction of the combined 
mechanism is 0.3 at /rep = 60 compared with a continuity frac- 
tion of zero for the cell -cell repulsion mechanism (figure 5c, 
pvalue = 1.1 X 10"^). 

The robustness of the combined mechanism to inhibition 
of the pinocytosis rate and to inhibition of the ceU-cell repul- 
sion strength disagrees with the experiments discussed above, 
in which inhibition of either mechanism reduced lumen 
formation. This discrepancy between our model and the exper- 
imental observations could be caused by the fact that the 
mechanisms of vacuolation and cell -cell repulsion are likely 
to be intertwined, while they are completely separated in the 
model. First, blockage of integrins to inhibit pinocytosis could 
also affect cell surface polarization, which is a crucial step for 
cell- cell repulsion. Second, cleavage or neutralization of the 
negative extracellular charged domains of CD34-sialomucins 
can affect the cytoskeleton, and thereby might also affect the 
transport and exocytosis of vesicles and vacuoles at the apical 
membrane. The intracellular domains of CD34-sialomucins 
are linked to the cytoskeleton, which disengages upon neutral- 
ization of the extracellular domains [45]. Thus, the discrepancy 
between the combined mechanism and experiments does not 
necessarily mean that the combined mechanism is incorrect, 
but could reflect the high level of cross talk between the 
mechanisms that drive vacuolation and cell -cell repulsion. 



2.9. Synergy of vacuolation and cell -cell repulsion 

We next asked whether the robustness of the combined mech- 
anism is due to synergy of vacuolation and cell -cell repulsion. 
Figure 5d shows the combined effect of the pinocytosis rate Ppi^ 
and the repulsion strength Jj-ep on the continuity fraction of 
the combined mechanism. Lumens are practically never con- 
tinuous for low values of Ppi^ and /rep (quadrant III), and 
practically always continuous for high values of both (quad- 
rant I). The high continuity fractions in quadrant IV indicate 
that vacuolation (dependent on Ppi^) reinforces cell -cell repul- 
sion (dependent on J^ep) in lumen formation. Similarly, the 
high continuity fractions in quadrant II indicate that cell -cell 
repulsion reinforces vacuolation in lumen formation. 

How does this synergy arise in the model? The cell -cell 
repulsion mechanism assists the vacuolation mechanism by 
breaking barriers, formed by adhering cells, between the sub- 
lumens. Additionally, cell -cell repulsion enlarges and 
stabilizes small sublumens that are created by secretion of 
vacuoles between cells. The vacuolation mechanism also 
reinforces the cell -cell repulsion mechanism. In the absence 
of vacuolation and for low cell -cell repulsion strength, cells 
do not detach, and the vessel thus remains solid. An increase 
in the pinocytosis rate (quadrant IV) results in the formation of 
sublumens, which drive cells to reposition into a multicellular, 
overlapping configuration. Cells that become overlapping can 
now detach by cell -cell repulsion in combination with 
secretion of vacuoles between the overlapping cells. Single 
cells that span the vessel can become pierced by vacuolation 
to create a continuous lumen. However, cell -cell repulsion 
for high repulsion strengths (figure 5d, /rep > 140) is more 
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Figure 5. Robustness and synergy, [a] Lumen formation by vacuolation at MCS 10 000 for different values of the pinocytosis rate {P^ J. [b] Continuity fraction of the 
combined mechanism and the vacuolation mechanism as a function of Ppjn. (c) Continuity fraction of the combined mechanism and the cell -cell repulsion mech- 
anism as a function of the repulsion strength (7rep). id) The continuity fraction of the combined mechanism expressed in a heat map as a function of Ppjn and ^p- (^) 
Lumen formation by cell-cell repulsion at MCS 10 000 for different values of V^ep- 



efficient in the absence of pinocytosis than for low pinocyto- 
sis rates. For the latter, cells do not get the chance to overlap, 
as will be explained in §2.10. 

2.10. Cells in unicellular sprouts need to reposition 
for cell -cell repulsion 

A priori, we had expected that cell -cell repulsion could only 
generate a continuous lumen in multicellular tubes and not 
in tubes with linearly, head-to-tail arranged cells such as in 
the tips of the branch. However, the cells in the branch tips 
repositioned to an overlapping (brick-like), multicellular con- 
figuration, in line with experimental observations [19,46], and 
this overlap allows lumen formation by cell -cell repulsion. 
We found that this repositioning is driven by two counter- 
acting forces between cells: strong attachment and active 
cell -cell repulsion. Cells do not immediately detach from 
each other upon active cell -cell repulsion as the standard 
CPM model treats the attachment of adjacent cell membranes 
as a vacuum. The membranes must zipper apart to let the 
fluid seep in. In the presence of vacuolation (figure 5d, 
/rep > 140), fluid is created between cells, and cells use 
cell -cell repulsion to instantly detach rather than to overlap. 

To check whether this zippering affects the efficiency 
of lumen formation by the cell -cell repulsion mechanism, 
we adapted the CPM to allow de novo insertion of ECM 



fluid [47]. With probability Pj, we consider the change in 
effective energy (AE) resulting from ECM fluid insertion 
(&CM medium (^) ^ ^1 ^ rather than from the extension of the 
cell membrane (4eii(^) ^ ^')- t^^is extended CPM, the orig- 
inal CPM is recovered for Pj = 0. Higher values of Pi allow 
cells to insert ECM fluid in between the repulsive apical mem- 
branes. Figure 6a, h shows magnifications of the simulations at 
MCS 1000 for Pi = 0 (figure 4a and the electronic supplemen- 
tary material, video SI) and for Pi = 0.1 (see electronic 
supplementary material, video S4), respectively. These magni- 
fications show that cells overlap in the tips of the branch for 
Pi = 0 and not for Pi = 0.1. Fluid insertion between cells thus 
seems to reduce overlap of cells, and thereby prevents lumen 
formation in the tips of the branch (see the electronic 
supplementary material, figure S2, for quantified results). 

We next studied the robustness of each mechanism to 
changes in the value of Pi. Figure 6c shows that the combined 
mechanism (green curve) is more robust to changes in Pi than 
the cell -cell repulsion mechanism (blue curve). The vacuolation 
mechanism (red curve) was not continuous for reference 
settings and this does not significantly change for higher 
values of Pi. Notably, the continuity fraction of the combined 
mechanism is often higher than the continuity fraction of the 
cell- cell repulsion mechanism and the vacuolation mechanism 
combined (e.g. for Pi = 0.3, pvalue = 1.1 x 10~^), another 
indication of reinforcement. 





Figure 6. Effect of cell overlap, {a) Enlarged view of overlapping cells for P\ 
fraction of all three mechanisms as a function of P\. 

2.11. Vacuolation requires impermeable vessel walls 

The vacuolation mechanism and the cell -cell repulsion mech- 
anism assume a different permeability of cells for fluid. Kamei 
et ah [3] showed that red quantum dots in lumens, formed by 
vacuolation, do not mix with the extracellular environment. 
The red quantum dots serially transferred from the dorsal 
aorta to previously unlabelled vacuolar compartments of 
ISVs in zebrafish. By contrast, the cell -cell repulsion mechan- 
ism assumes open connections through paracellular openings 
and fluid can flow into the lumen from the ECM [17]. To 
mimic the effect of permeable vessel walls, we tested the behav- 
iour of the model for Afi^ids = 0. In this case, the fluids become 
completely 'compressible'. Only the vacuolation mechanism 
requires fluid incompressibility for generating continuous 
lumens (see electronic supplementary material, figure S3), 
because the lumens otherwise continuously coUapse. Thus, our 
simulations suggest that vessel walls must be impermeable to 
water (or lumen fluid must be actively replenished to maintain 
hydrostatic pressure) in order to generate continuous lumens 
by the vacuolation mechanism, but this is not required for the 
cell -cell repulsion mechanism or for the combined mechanism. 



3, Discussion 

Extensive experimental research has resulted in two alternative 
proposed mechanisms of lumen formation: vacuolation [3,4] 
and cell-cell repulsion [5,17]. Our computational model 
suggests that the two mechanisms may act synergistically in 
lumen formation. Cell -cell repulsion can reinforce vacuolation 
by stabilizing sublumens and by separating cells to connect 
sublumens with the ECM fluid. Vacuolation can reinforce 
cell -cell repulsion by creating sublumens, which repositions 
cells into an overlapping configuration, and by piercing 
single cells in the tips of the branch. Additionally, vacuolation 
could assist in the expansion of the lumen, which was so far 



= 0. [b] Enlarged view of reduced overlapping capacity for P\ = 0.1. (c) Continuity 



suggested in the cell -cell repulsion mechanism to occur by 
cell shape changes. 

We validated our model assumptions and simulation results 
based on published experimental evidence. As previously dis- 
cussed, for lack of quantitative values of the model parameters, 
the model can only make qualitative predictions. First, without 
apical -basolateral cell surface polarization, no lumens are 
formed in our model (figure 2h). This agrees with experiments 
in which polarization was prohibited in the absence of functional 
vascular endothelial (VE)-cadherin or phosphatase and tensin 
homologue [6,17]. Second, in our simulations of the vacuolation 
mechanism, lumens cannot form for low pinocytosis rates 
(figure 5h). Experimentally, lumen formation is indeed pre- 
vented by blockage of integrin signalling for pinocytosis 
[18,44]. Third, in our simulations of the cell-ceU repulsion mech- 
anism, continuous lumens cannot form at low repulsion 
strengths (figure 5c). In agreement with this model prediction, 
neutralization or cleavage of the extracellular negative domains 
of CD34-sialomucins reduces lumen formation [5]. 

Although the model simulations suggest that the vacuo- 
lation and repulsion mechanisms act synergistically, one may 
ask whether both mechanisms indeed co-occur in vivo. It is 
possible that vacuolation and cell -cell repulsion function 
in different types of vessels; vacuolation for single-cell 
capillaries (e.g. ISV of zebrafish) and cell-cell repulsion for 
multicellular tubes (e.g. dorsal aortae of mice) [3,6,9-12]. There- 
fore, we tested how the vacuolation, cell -cell repulsion and 
combined mechanisms would act in different t3rpes of vessels, 
which we represent in our model by initial configurations of 
one, two and three layers of cells (figure 7). Cell -cell repul- 
sion is not functional in a one-cell-thick vessel with aligned 
cells (figure 7a). Cells do start to overlap, but not sufficiently. 
Vacuolation does create a lumen, but it continuously collapses 
again (figure 7h). The combined mechanism forms stable 
lumens (figure 7c). In two-cell-thick vessels cell -cell repulsion 
is much more efficient than vacuolation (figure 7d-f) and 
stable lumens form by cell-cell repulsion and by the combined 
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Figure 7. Lumen formation by each mechanism in different vessel types, (a-i) Lumen formation by cell -cell repulsion, vacuolation or the combined mechanism 
from left to right for an initial vessel thickness of one cell (top row), two cells (second row) or three cells (third row) at MCS 10 000. (j,k). Time-lapse images (MCS 
500, 1000 and 2000) of lumen formation by the combined mechanism in a one-cell-thick vessel and in a two-cell-thick vessel. 
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mechanism, but not by vacuolation. In three-cell-thick vessels, 
the cell -cell repulsion and combined mechanisms can repro- 
duce cavitation, the apoptosis of cells in the middle of the 
vessel that detached from the vessel w^all (figure 7g,z), w^hich 
is often seen in epithelium [13]. Interestingly, the combined 
mechanism produces phenomena that optically resemble 
aspects of vacuolation or of cell -cell repulsion, depending on 
w^hether it acts in one-cell-thick vessels or in tw^o-cell-thick 
vessels. The optical resemblance to vacuolation is illustrated 
in figure 7/ and the resemblance to cell -cell repulsion in 
figure 7k, w^hich show^ the temporal development of lumen for- 
mation by the combined mechanisms in a one-cell-thick vessel 
and a two-cell-thick vessel, respectively. Thus, although vacuo- 
lation is observed in capillaries and cell -cell repulsion in 
multicellular tubes, the combined mechanism could be the 
underlying mechanism in both vessel t3rpes. 

To validate w^hether the combined mechanism is actually 
functional in vivo, we propose the follow^ing experiments. 
Based on our modelling results in figure 7, we expect that 
inhibition of vacuolation by reducing pinocytosis should pre- 
vent lumen formation in capillaries (as in figure 7a), but not 
in larger vessels (as in figure 7d,g), and inhibition of cell- 
cell repulsion by cleavage of negatively charged extracellular 
proteins should destabilize lumen formation in all vessel 
types (as in figure 7h,e,h). Here, simply visualizing vacuoles 
microscopically w^ill not suffice for validating the combined 
mechanism: lumens also form w^ith cell -cell repulsion in 
combination w^ith secretion of pinocytotic vesicles, w^ithout 
fusion of vesicles into vacuoles, depending on the amount 
of vesicles (see electronic supplementary material, figure S4). 

Besides synergy of the vacuolation and cell -cell repulsion 
mechanisms, our model also predicts that lumen formation 
by cell -cell repulsion in unicellular tubes requires cells to 
reposition into a (brick-like) overlapping, multicellular con- 
figuration. Our model results suggest that cell -cell adhesion 
along w^ith strong cell -cell repulsion strengths facilitate the for- 
mation of such overlapping cellular configurations. To validate 
this mechanism, it should be established experimentally if cells 
immediately detach once CD34-sialomucins are in the apical 
membrane, or that strong adhesion keeps them attached for 
time spans sufficiently long to allow repositioning of cells. 

These validation experiments all result from qualita- 
tive model predictions, as the quantitative values for most 



parameters are not know^n. Dose-dependent experiments 
for cell adhesion strengths, pinocytosis rates and cell -cell 
repulsion strengths can help to tune the representing quali- 
tative parameters in the model. Additionally, various 
parameters could be quantified experimentally to allows for 
quantitative model predictions. The adhesion strength of 
cells can be quantified by the force that is required to pull 
them apart [48]. This method might also be useful to find 
quantitative values, or at least the relative ordering, of the 
contact energy parameters that describe adhesions bet- 
ween polarized cells, non-polarized cells and possibly also 
the ECM. Other contact energy parameters of our model 
describe adhesions of subcellular compartments such as ves- 
icles, w^hich cannot be quantified in this way. Instead, vesicles 
and vacuoles can be visualized microscopically and their 
speed and type of movement (e.g. diffusive) could thus 
be quantified. 

Our model can become a useful tool for designing 
new experiments and new insights into lumen formation. 
We propose three new research questions in which we 
believe cooperative computational and experimental research 
is important. First, what is the exact function of several key 
proteins in lumen formation? Cdc42 and Moesinl are, for 
instance, suggested to be involved in polarization of the 
cell, in targeting of vesicles to the apical membrane and in 
structural changes of the cytoskeleton for cell shape changes 
[6,17,29]. It is difficult to pinpoint their exact function by 
experiments only, because lumen formation fails all together 
in the absence of these proteins. Second, how is lumen for- 
mation regulated in dynamically growing sprouts? To focus 
on lumen formation, we started our model with a preformed 
sprout. To gain insights into the regulation of angiogenesis 
as a whole, the model could be extended with ECM re- 
modelling and dynamic sprouting. Third, if the combined 
mechanism, indeed, drives lumen formation, then how are 
the two mechanisms regulated and balanced to locally opti- 
mize lumen formation? For each question, the model can 
be used to test consistency of hypotheses, which can pro- 
vide new insights and help to guide new experiments. 
In conclusion, in collaboration with experimentalists, our 
simulation model can contribute to a better understanding 
of the mechanisms of lumen formation during blood vessel 
development. 
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Figure 8. Representation of cells in the compartmental CPM. The compart- 
mental CPM projects cells on a regular lattice and cells can consist of multiple 
compartments. All lattice site, x, of the same cell are outlined by a bold black 
line, representing that they share the same unique cell identifier ^{cr{x)). 
Each lattice site has a number that represents its compartment identifier 
cr{x) and its type T{cr{x)) is indicated by its colour. 



4. Material and methods 

We developed an agent-based, computational model of lumen for- 
mation that connects the subcellular, cellular and the Vessel' 
scales. The CPM [24,25] describes the motility, shape and physical 
interactions of cells. To model polarization of the cell surface and 
the creation of fluid-filled vesicles and vacuoles, we use an exten- 
sion of the CPM in which cells can compartmentalize [26-28]. In 
this section, we first explain the extended version of the CPM. 
Next, we outline the modelling of pinocytosis, vesicle and vacuole 
movement and secretion. 



4.1. Compartmental cellular Potts model 

The CPM projects cells on a regular lattice (figure 8). Each lattice 
site, X, is associated with a unique compartment identifier a{x) 
and has a type t((t(x)), which can be cytoplasm, apical, hasolateral, 
vesicle, vacuole, ECM, ECM fluid or luminal fluid (figure 3). Initially 
cells, ECM and ECM fluid consist of a single compartment with 
a unique cell identifier f(cr(x)). Additional compartments, which 
are formed in a cell upon membrane polarization and vacuo- 
lation, obtain the same cell identifier f(cr(x)). New cell 
identifiers with type luminal fluid are created upon secretion of 
vesicles and vacuoles. 

To model random motility, fluids and subcellular compart- 
ments (except for vesicles) move by copies from x to a randomly 
selected second neighbouring lattice site x' (figure 3a). A copy 
changes the effective energy of the system (E), which depends on 
the areas and contacts (adhesive properties) of cells, compartments 
and fluids: E = Earea + Econtact- A copy is always accepted for 
a decrease in effective energy and is otherwise accepted with a 
Boltzmann probability PBoitzmann(E) = e^~^^^^\ where ijl resembles 
active cell motility. 

Earea IS the costs iu energy for deviation of the actual area (a) 
from the preferred target area (A), with an elasticity lambda, and 
is given by 

Earea = Aarea(f) (A(f) - a{^)f + Aarea(o-) (A(o-) - a{cT)f. 

We model fluids as near-incompressible fluids; fluids are cell 
identifiers of type ECM fluid or luminal fluid with Aarea( f ) ^ 0 

and X^reaH / 0. 

Contact energy (Econtact) mimics adhesion and repulsion at 
interfaces of compartments. There are two types of contact 
energy: internal (Ji) and external (/e). Internal contact energy is 
defined between compartments of the same cell, and external 



contact energy between compartments of different cells. The 
total contact energy is defined as 

£co„tact = E(..^)/E(^('^W).^('^(*')))(l-5(a<^W),a<^(^)))) 

with the Kronecker delta 8(x, y) = {1, x = y; 0, x / y}. All 
/-values are listed in tables 1 and 2. Cells (non-polarized or polar- 
ized) stay connected in the absence of lumen formation 
mechanisms (figure 2h,c) for reference /-values. 

4.2. Pinocytosis 

During each MCS, as many copy attempts as there are pixels in the 
lattice are performed (550 x 550). Extra mechanisms (m) of lumen 
formation are performed after every n^ MCS (figure 9). Cell surfaces 
polarize every other MCS (figure 9; mechanism 5; ns = 2) to aUow 
these polarized membrane pixels to internalize by usual CPM 
copies, representing invagination of ECM fluid at the membrane 
as seen during pinocytosis. Such internalized polarized membrane 
pixels become compartments of type vesicle with probability Ppi^ or 
otherwise part of the cytoplasm (figure 9; mechanism 6; = 2). A 
vesicle is kept one pixel in size by a target area of one and a high 
lambda. As ECM fluid is taken up into the vesicle during pinocyto- 
sis, the target area of ECM fluid decreases, while the target area of 
the pinocytosing cell increases by one. 

4.3. Vesicle and vacuole movement 

Vesicles swap position with a randomly selected neighbour once 
per MCS (figure 9; mechanism 1; ni = 1). This swapping is per- 
formed with a pre-set acceptance probability Pa multiplied by 
^'Boitzmann(E), with E the resultiug effective energy of the swap. 
Pa tunes the velocity of the vesicle. 

Fusion events of vesicles and vacuoles can occur every MCS 
(figure 9; mechanism 2; n2 = 1) and during a copy of a vacuole 
over a vesicle. Every MCS, neighbouring compartments of type 
vesicle or vacuole fuse with probability Pfuse ir^to a single compart- 
ment with type vacuole. The target area of the formed vacuole is 
the sum of the target areas of the fused compartments. Similarly, 
when a vacuole copies over a vesicle of the same cell, the target 
area of the vacuole is increased by the target area of the vesicle. 

As vesicles can become vacuoles by fusion, small vacuoles 
can also become vesicles. First, single pixels of type vacuole that 
are split off from a vacuole, called the donor vacuole, and are 
surrounded by cytoplasm become vesicles each MCS (figure 9; 
mechanism 3; n^ = l). Second, a donor vacuole that became a 
single pixel in size by a copy becomes a vesicle. In both cases, 
the created vesicle gets a target area of one, unless the donor 
vacuole in question had a target area of zero. Then, the vesicle 
is assigned with a target area of zero and will soon be deleted 
by regular CPM movements. If a donor vacuole remains, then 
the target area of the created vesicle is subtracted from the 
target area of the donor vacuole. If a donor vacuole becomes 
one pixel in size by a copy event and has a target area larger 
than one, then its residual target area must be redistributed to 
remain a target area of one for the created vesicle. If the donor 
vacuole was copied over by luminal fluid or a vacuole, then 
the residual target area is added to that compartment. Other- 
wise, the residual target area is added to ECM fluid to keep 
the total target area of the system constant. 

4.4. Secretion 

Secretion can occur every other MCS (figure 9; mechanism 7; 
n7 = 2) and by a copy of a fluid over a vesicle or vacuole. 
Every other MCS, a pixel at the membrane of type vesicle or vacu- 
ole, together with all first-order connected pixels of type vesicle or 
vacuole, becomes a compartment of type luminal fluid. The 
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Figure 9. Flow chart of the model. During each time step, also called Monte Carlo step (MCS), as many copy attempts as there are pixels in the lattice (550 x 550) 
are attempted. After each time step, or after every other time step, additional functions for lumen formation are performed. 



combined target area of the simultaneously secreted vesicles and 
vacuoles is assigned to the luminal compartment and is sub- 
tracted from target area of the secreting cell. If only a part of 
the vacuole is secreted, then the size of the secreted part is sub- 
tracted from the target area of the vacuole, leaving a vacuole with 
a minimal target area of zero. Secretion can also occur by a copy 
X to x' , for v\^hich T{a(x')) is vesicle and T{a(x)) is luminal fluid 
or ECM fluid, resulting in a decrease of one of the target area of 
the secreting cell and an increase of one of the fluid. Similarly, the 
target area of the secreting cell decreases by one w^hen its vesicle 
is copied over by a compartment of the same cell of type apical, 
basolateral or cytoplasm or by a compartment of another cell. To 
conserve the total target area of the system, the target area of 
ECM fluid is then increased by one. 

ECM fluid and luminal fluids can fuse every MCS (figure 9; 
mechanism 4; n4 = 1) and by some copy events. Neighbouring 



cells of type luminal fluid are fused every MCS and can fuse to 
ECM fluid wrhen the luminal fluid is in contact w^ith the sur- 
rounding ECM or w^hen it is completely surrounded by ECM 
fluid. When luminal fluid copies over luminal fluid of a different 
cell identifier and thereby deletes it, the target area of the latter is 
added to the first. If luminal fluid is deleted by a cell type other 
than luminal fluid, then its target area is added to ECM fluid. 
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